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Abstract 

Time-dependent Hamiltonian dynamics is derived for a curve (molecular strand) in R 3 that experiences both 
nonlocal (for example, electrostatic) and elastic interactions. The dynamical equations in the symmetry-reduced 
variables are written on the dual of the semidirect-product Lie algebra so(3) (S) (R 3 © R 3 © R 3 © R 3 ) with three 
2-cocycles. We also demonstrate that the nonlocal interaction produces an interesting new term deriving from 
the coadjoint action of the Lie group 50(3) on its Lie algebra so(3). The new filament equations are written in 
conservative form by using the corresponding coadjoint actions. 

To cite this article: D. D. Holm and V. Putkaradze, C. R. Acad. Sci. Paris, Ser. I XXX (2008). 

Resume 

Nous derivons la dynamique Hamiltonienne d'une courbe (chaine moleculaire) dans l'espace physique R 3 sujette 
a des interactions elastiques ainsi que non- locales (electrostatiques par exemple). Les equations dynamiques des 
variables reduites par symetrie sont ecrites sur l'espace dual de l'algebre de Lie so(3) (S) (R 3 © R 3 © R 3 © R 3 ) 
(produit semidirect) avec trois 2-cocycles. Nous demontrons aussi que l'interaction non-locale produit un nouvel 
terme interessant, qui derive de Taction coadjointe du group de Lie 5*0(3) sur son algebre so(3). Les nouvelles 
equations du filament sont ecrites sous une forme conservative grace aux actions coadjointes correspondantes. 
Pour citer cet article : D. D. Holm and V. Putkaradze, C. R. Acad. Sci. Paris, Ser. I XXX (2008). 



Version francaise abregee Cette note depasse le cadre des approches a la Kirchoff pour deriver la dy- 
namique Hamiltonienne d'une courbe (chaine moleculaire) dans l'espace physique R 3 , lorsqu'elle est 
sujette a des interactions elastiques ainsi que non- locales (electrostatiques par exemple). Cette note 
s'inspire d'une extension de la theorie des barres geometriques parfaites [12], qui presente l'interaction 
elastique dans un contexte geometrique par l'usage des coordonnees materielles. Cette theorie a deja 
ete utilisee pour la description de quelques aspects de la dynamique des proteines [2]. Toutefois la ge- 
neralisation des ces theories aux interactions non-locales (qui dependent du temps) exige l'application 
de methodes geometriques pour deriver le principe d'action (15), qui vont au-dela de l'approche a 
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la Kirchoff. Les equations dynamiques des variables reduites par symetrie (16,17) sont ecrites dans 
l'espace dual de Palgebre de Lie so(3) (S) (M 3 ©]R 3 ©R 3 ©IR 3 ) (produit semidirect) avec trois 2-cocycles. 
De plus, les nouvelles equations du filament sont ecrites sous une forme conservative (25) en utilisant 
les actions coadjointes correspondantes. Notre approche rend possible l'inclusion coherente des effets 
electrostatiques et inertiels dans les etudes theoriques et numeriques de la dynamique des chaines 
biologiques. Cela garantit aussi la conservation de l'energie et etend les traitements Hamiltoniens a 
la Lie Poisson [12,6,10] afin d'englober la dependance non-locale des variables. 

1. Introduction 

In contemporary science, two basic approaches are taken in describing the dynamics of biological 
strands (such as proteins). These two approaches might be termed molecular dynamics and geo- 
metric dynamics. Molecular dynamics treats a biological molecule as a collection of charged masses 
interacting in force fields; the equations of motion are derived from Newton's Second Law, with 
proper addition of random forces. Using this method, scientists have been able to accurately model 
the realistic dynamics of proteins with complex shapes [13]. Even though this approach has been 
successful, it poses tremendous demands on computing power; it also limits the ability to achieve 
theoretical understanding of the dynamics. In addition, the computation of folding of even the most 
basic proteins using such detailed dynamics challenges the limits of most modern computers. This is 
due to the necessity of modeling the motion of every single atom in the protein as well as using the 
atomistic stochastic approach for computing resistance and self-interaction of the molecular strand 
in water, or other solvent. 

Alternatively, consideration of a protein or biological strand as a continuum curve has a long and 
rich history with some notable successes [8,1,7,9]. Most considerations in this approach - both in the 
investigation of stationary and time dependent models - have been made in the framework of the 
Kirchhoff model of elastic rods. That framework, however, has been so far not been able to overcome 
the mathematical difficulties of incorporating electrostatic (nonlocal) effects into the dynamics, by a 
straight-forward extension of Kirchhoff theory. Some recent progress has been made on the stationary 
solutions in the continuum framework [5]. However, no time-dependent theory considering nonlocal 
effects yet exists to our knowledge. The difficulties arise because Kirchhoff's theory is formulated 
in the intrinsic frame connected with the deformed rod. Hence, the computation of distances in 
Euclidean space becomes non-trivial, thereby increasing the difficulty of constructing a consistent 
time-dependent theory. 

This paper overcomes the problems confronted in standard Kirchhoff-based approaches by using an 
extension of the theory of exact geometric rods [12], which puts elastic interactions into a geometric 
framework using the material frame viewpoint. This theory has already been used to describe some 
aspects of protein dynamics [2]. However, the generalization of these theories for time- dependent 
nonlocal interactions requires geometric methods that go beyond the Kirchhoff approach. The paper 
considers only inertial effects. The introduction of dissipation into the motion of biological strands 
is a complicated issue [8,7] that will be deferred to future work. 

Acknowledgements. The authors were partially supported by NSF grant NSF-DMS-05377891, the 
US Department of Energy, Office of Science, Applied Mathematical Research, the Royal Society of 
London Wolfson Research Award and the MISGAM program of the European Science Foundation. 

2. Motion of exact self-interacting geometric rods 

We consider rigid conformations of charges mounted along a flexible filament at r(s,t) at distances 
rj k (s, t) and allow these charges to interact with each other via a potential (for proteins, the screened 
electrostatic potential). For simplicity, each charge 'bouquet' is assumed to rotate as a rigid body with 
respect to its origin. This rigid conformational rotation is illustrated in Figure 1. The dependence 
on time t of the orientation of the rigid conformation of charges at a spatial point r(s, t) along the 
filament is denoted as A(s,i) G SO(3). Suppose each rigid conformation of charges is identical and 
the k-th electrical charge is positioned near a given point r(s,t) of the curve at the reference state 
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Figure 1. Rigid conformations of charges are distributed along a curve. 

r+r] k (s). Here, r) k (s) is a vector of constant length that determines the position of the A:-th electrical 
charge relative to the point r(s, t) along the curve in its reference configuration. If the curve position 
r remains fixed, rotation is allowed only in a plane so A G 50(2) and there is only one charge k = 1, 
our model reduces to that considered in [11]. In general, the position c k of the fc-th charge in the 
rigid conformation anchored at position r(s,t) rotates to a new position with A(s, t) £ 50(3) as 

Cfc(s,t) = r(s,t) + A(s,t)r) k (s) , where A(s,0)=ld. (1) 

The key to further progress is to use maximum possible reduction of the Lagrangian to the 50(3)- 
invariant quantities. As we show below, because of the nonlocality, the complete reduction from the 
Lie group to its Lie algebra is impossible, yet the equations can still be formulated as motion on a 
Lie algebra with an elegant mapping from the unreduced Lie group terms to the Lie algebra. 
The nonlocal part of the potential energy l np of interaction between rigid conformations of charges 
at spatial coordinates s and s' along the filament depends only on the distance dfc, m (s, s') between 
the fc-th and m-th charges in the two conformations, 

l n p = ^2^ J u(dk, n (s,s')jdsds' , dk, m (s,s') = \c k (s) -c m (s')\ . (2) 

fc,m 

This scalar distance may also be expressed in terms of vectors seen from the frame of orientation of 
the rigid body at a point x along the filament, as 

dk, m (s,s') = |c fc (s) -c m (s')\ = \k(s,s') +r) k (s) - f(s, s')r] m (s')\ , (3) 

where k(s, s') = A~ 1 (s)(r(s) - r(s')) E R 3 and £(s, s') = \- 1 (s)A(s') £ SO (3) . (4) 

The first of these quantities is the spatial vector r(s) — r(s') between two points on the filament, 
as seen from the orientation A(s) of the rigid body (charge conformation) at s on the filament. 
The second is the relative orientation of the rigid bodies (charge conformations) at s and s' . A 
transposition identity, £(s, s') T = £(s' , s) = £(s, s') _1 j follows from the definition of £(s, s') in (4). 
Remark 1 Both the vector k(s, s') and the relative orientation £(s, s') defined in (4) are invariant 
under changes of the orientation of the spatial coordinate system obtained by the left action r(s) — 
r(s') — > 0(r(s) — r(s')) and A — > OA of any element O in the rotation group 50(3). 
The localized part of the Lagrangian depends on A, A, A', r(s), r, r' , where the dot denotes the 
time derivative and prime is the derivative with respect to arclength s. If it is left-invariant under 
the action of 50(3), this Lagrangian may be reduced to a function of left-invariant quantities: 
Q = A _1 A' £ so(3), w = A- 1 A £ so(3), T = A~V E R 3 , 7 = A" 1 /- e R 3 and p = A _1 r £ R 3 . 
Capital Greek letters denote derivatives in s, while lower-case Greek letters (except for p) denote 
derivatives in time. Bold letters such as u> denote vectors in R 3 while CI is a skew 3x3 matrix in 
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the Lie algebra so(3) whose entries correspond to vector components via the isomorphism between 
so(3) and R 3 . For example, for any vector u£l 3 one has flv = x v. In components, this is the 
map Ctjk = — tjki^i without any extra adornments. 

In terms of these quantities, we assume the symmetry-reduced Lagrangian L may be written as the 
sum of a local part I and a nonlocal part l np , according to 

L(p,7,r,w,n,o = i(p, 7, r,u> > n)+ / u( K (s,s'),fa,s'))dsds' -.= i + i np , (5) 

where the left-invariant quantity rc(s, s') in (4) may be expressed using p and £ as 

k(s, s>) := A- 1 ( S ) (r(s) r(s')) = p(s) fa, s')p(s') . (6) 

Kinematics. Let us compute the space and time derivatives of p = A _1 r € R 3 . The space (s) 
derivative of p (denoted by a prime) and time (t) derivative (denoted by a dot) are given by 

p' = — Op + r = -f2xp + r and p=-wp + j= -ujxp + j. (7) 

Compatibility of these formulas arises from equality of the cross-derivatives of r and A. Namely, 

f +wxr=7' + fix7 and tl = fl X u; + u; ' . (8) 

Variations. The variations of p, u>, 7, fi and T are computed by the following steps: 

5p = -A^SAA^r + A _1 <5r = £p + * = -£xp + * = px£ + *, (9) 

where one defines the left-invariant variations £ = A _1 <5A G so(3) and \I/ = A _1 <5r £ R 3 . The 
variations £7, 8T, 5ui and 8Q are given in terms of the the left-invariant quantities £ and \& by 

a* a* 

57 = - X x 7 + w x * + — = 7 x S + u> x * + — , (10) 

ff = -sxr + fix* + — =rxE + fixfi- , (11) 

os os 

5u = uxS + — and = fl x £ + — . (12) 
at os 

The key to understanding the nonlocal variations lies in the matrix formula 

tT x 8fa,s') = -Ad e -i (SjS /)S(s)-|-S(s'), where Ad s -iS := £ _1 X;£ , (13) 

obtained from the definition of £(s, s') in equation (4). The variation of k in (6) is then given by 

8k(s, s') = - £(s) x k(s, s') + *(s) - f (s, »')*(*') ■ ( 14 ) 

Calculation of energy variations. The equations of motion are computed from the stationary 
action principle SS = 0, with S = J Ldt and L = I + l np in equation (5), for which 

*, 1 „) + ^,*.) + ^ r ^, r . 4e ) <B _0. (15) 

The terms proportional to £ and \P give, respectively, 

9 5/ 5l\ f d 61 ^ 61 \ 51 61 81 

+ u; x — + h n x — = — X7H xTH xp (16) 



J ^{s,s')xK(s,s')ds , + J Z{s,s')ds', 
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d 61 81 \ ( d 51 n 51 

di^ + UJX 5^) + {d- s 5f + nX 5f 



- + j- M -t M -v,swy (it) 



Here, we have defined the nonlocal contribution 

,ATT X T fffj 

Z(s,s')~Z( S ,s')x=Z(s,s')(^r( S , S 'y) -|i( s , s ')£ T ( S , S '). (18) 

The term Z(s, s') is the contribution from the nonlocal part of the Lagrangian we have sought. A 
direct calculation shows that Z T = —Z, so Z G so(3). This expression appears naturally in geometric 
mechanics of interacting oriented bodies like asteroids [3] , and more general theory of reduction by 
stages [4] illuminates the geometric structure of such interactions. 

3. Spatial conservation laws 

The physical meaning of equations (16, 17) is revealed by writing them as spatial conservation laws. 
For this, we invoke the following identities valid for the Ad- and Ad*-actions of any element g(t) € G 
in the Lie group on r\ £ g in the Lie algebra and /i £ g* in its dual with pairing ( • , • ) : g* x g — ► K 

Mi Ad 9- I (*)^ Ad ff(t) ? ?) = (j"> ad ^ 7 ?) aild ( Ad S« j^ Ad 9- 1 (*)' U ' 7 ?) = ( ~ ad ^' 7 ?) ' ( 19 ) 

where Ad* : G x g* — > g* is defined by (Ad*/ t \/U, 77) := (/i , Ad fl ( t )7j) and a = g g'(t) 6 g belongs 
to the Lie algebra g. Equations (16, 17) are formulated on the dual of the Lie algebra, for which the 
second equation in (19) gives 

Ad; (t) ^(Ad;_ 1(4) /x(t)) , ri) = ^i-adUhv)- (20) 

To derive the conservation form of equations (16, 17), we consider the group G = ££7(3) with the 
group element g = (A(s,t),r(s,t)), whose left-invariant time-derivative is a = (A _1 A,A _1 r) = 
(u>,7). Using the definition of ad* for the se(3) Lie algebra yields 
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» ,* f til til \i d f 51 51 \ / 51 51 51 \ 

Ad^^Ad* I(t) — ,— =- _ — + wx — + 7 X— ,-wx— . (21) 



The tangent vector obtained from the derivative in arclength s is (A~ 1 A / , A^V) = (0,r). Thus, 



Ad 



f 51 5l\-\ Of 01 01 \ / 01 01 _ 01 \ 

I . — = — I . — U Ox kTx — -Ox — 



9_1 WVjn ' ST J 



d t til 51 \ / 61 51 „ 61 



.Ox — + TX— ,-Ox— . (22) 
dsVjO' 5VJ V 50 6T' SrJ y J 



Remark 2 The nonlocal term (18) arises as the derivative of the nonlocal part of the potential with 
respect to Lie algebra elements ft and T, as follows. 

Upon identifying coefficients of the free variations Ex = A _1 <5A and *Sf = A~ 1 5r, one may write the 
following identity relating different variational derivatives of the nonlocal potential l np : 



Using expressions (13) for (14) for 5k, (12) for 6ft and (10) for 6T, then collecting terms 

proportional to the free variation £ yields the following identity, which implicitly defines 5l np /6fl 
in terms of known quantities, 

d tilnp t ^ 51 np f dU l 
ds 6ft 

where we have defined Z(s, s') according to (18). Likewise, identifying terms multiplying gives 



d 5l np 
ds 5T 



nX lff = J S (S ' S ' ] X K(S ' S ' ] d3 ' + J Z(S ' S>) d3 ' ' {23) 
ccording to (18). Likewise, identifying terms multiplying gives 
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Therefore, we conclude that equations (16, 17) are equivalent to the following pair of conservation 
laws expressed using variations of the total Lagrangian L — I + l np : 



dt 



» j* fSL SL 



d_ 

ds 



Ad* 



/SL 5L\ 

-^Asn' sf) 



Ad! 



/SL SL\ 



(25) 



The terms (SL/Sp x p, 5L/5p) describe the influence of the external torques and forces, respectively, 
which are not expressible in conservation form. (This is similar to the situation for the heavy top.) 
Remark 3 On Legendre transforming the total Lagrangian L to the Hamiltonian, 



H(ji,0-,n,T, P ) = J (ii-u + p-i)da-L(u, T ,n,r,p) 

equations (7), (8), (16) and ( 17) may be expressed in Lie-Poisson form as 



(26) 



dt 







pX 


px 


(d s + nx) 


rx 


fix 




SHjSp = u) 


P 




px 











Id 




5H/Sp 


n 




(d s + flx) 
















SH/sn 


r 




rx 











(ds + nx) 




SH/ST 


P 




fix 


-Id 





(d s + nx) 







SH/Sfi = 7 



(27) 



This Lie-Poisson Hamiltonian matrix is dual to the semidirect-product Lie algebra so(3) (S) (R 3 ©K 3 © 
K 3 ©K 3 ) with three different types of 2-cocycles defined on its normal M 3 subalgebras. The symplectic 
2-cocycle in {p, (3} induces the generalized 2-cocycle in {T, f3}, for which ft is a Casimir; so fix is 
a constant in this 2-cocycle. In contrast, the quantity fix in the generalized 2-cocycle for {p, ft} is 
a connection form. The latter 2-cocycle also appears in the theory of complex fluids [12,6,10]. 
Remark 4 In this paper, the influence of nonlocality (e.g., charge screening) on rod mechanics was 
studied using the Euler-P oincare variational method. This variational approach led to an equivalent 
Lie-Poisson Hamiltonian formulation of the new equations (16, 17). Applying the Ad* g _-L^. trans- 
formation from body to spatial variables in these equations produced a great economy of form and 
exposed the meaning of the interplay among their various local and nonlocal terms. 
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